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Abstract 

In the framework of local thermodynamic equilibrium by Buchholz, 
Ojima and Roos, a class S x of observables, whose members are sup- 
posed to model idealized measurements of thermal properties of given 
states at spacetime points x, plays a crucial role in determining and 
characterizing local equilibrium states in quantum field theory. Here 
it will be shown how elements from this space can be reproduced by a 
specific model of the idealized measurements modeled by an Unruh-de 
Witt detector. 



1 Introduction 

In [4] a method to identify states in relativistic quantum field theory that 
allow locally a thermodynamic interpretation without being in global ther- 
modynamic equilibrium has been proposed. The basic procedure is the 
following: First one fixes at each spacetime point some (linear) space of ob- 
servables, denoted by S x , which model idealized measurements of thermal 
properties in arbitrarily small spacetime regions. Next one chooses a set of 
reference states with known thermal properties (in the cases investigated so 
far this have been mixtures of global equilibrium states). A given state is 
then called (locally) thermal at a spacetime point x if there exists a reference 
state such that the expectation values of all elements in S x in the two states 
agree, i.e. as far as local measurements of (expectation values of) thermal 
parameters are concerned, the given state looks like the reference state. 
If a state fulfills this condition of thermality at x one can then consistently 
assign to it the thermal parameters whose local measurement are modeled 
by elements of S x and among them the relations (e.g. equations of state) 
established for the reference states also hold true. As the reference state 
may vary from point to point, this yields for states fulfilling the thermality 
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condition in a region O an assignment of thermal parameters depending on 
x £ O (e.g. a time- and space-dependent temperature). For further details 
and applications of the formalism to models see [3], [1] and [5]. 
If the set of reference states is fixed (e.g. as mixtures of global equilibrium 
states), there remains the question of which tS^-spaces one should pick. In 
the models considered so far, the 5 x -spaces have been chosen as the spaces 
generated by so called "balanced derivatives" of Wick-squaresS , i.e. by 
elements 

lim d» (4>(x + CMx - C) - Woo {<t>{x + 0<t>{x - 0) i) 

(fi a multiindex, 1 the identity operator and lv^ the vacuum-state). Now a 
justification of this choice of S^-spaces is given in the models by the fact that 
it allows for a local assignment of the interesting thermodynamic properties 
to non-equilibrium-states in the above sense while still allowing non-trivial 
local-equilibrium states (i.e. states where the thermal parameters do vary in 
space and time). There is also a general discussion on the specification of a 
maximal space T x of "pointlike" observables that may be used to determine 
local thermal properties of states [4]; the above set 5^ should be a proper 
subset of T x . 

As the concept of local thermodynamic equilibrium states under discussion 
is based on the idea of idealized measurements at a point, one could hope 
that a more detailed description of such an idealized measurement in a model 
could give some additional insight into the choice of S x spaces. That this is 
indeed true will be shown in this article. 

As a model for the measurement process some kind of Unruh-de Witt detec- 
tor [TT],[7], in this case a two-level system moving along a given trajectory 
in spacetime and interacting with the given quantum-field, is used. For 
large interaction-times and time-invariant states the (suitably normalized) 
probability for the transition of the detector system from its initial (ground- 
or excited state) to its respective other state due to interactions with the 
quantum field can be used to determine thermal properties of the quantum 
field. Namely, by the principle of detailed balancing in a thermal state one 
expects the transition rates from the ground- to the excited and from the 
excited- to the ground-state to be related by a Boltzmann factor e~P E / h if 
the state of the field is a thermal state at inverse temperature (5 and E is 
the energy difference between the two levels of the detector system. In fact, 
such a relation between the two transition rates for all two-level monopole 

1 for curved spacetimes the "general covariant Wick-square" 8 , 2 can be used for a 
similar definition 



2 



detector "at rest" wrt. the thermal statcl exactly corresponds to the KMS 
condition for this state [TO]. More generally, the dependence of the transi- 
tion probabilities on E/H corresponds to "spectral properties" of the field, 
with the principle of detailed balancing as a relation between the rates for 
E and —E as a special case. 

Now in order to be able to proceed to measurements taking place in a short 
time-interval, the idea is first not to look at the absolute transition proba- 
bilities but rather at their differences to those in a common reference state 
(i.e. to "remove the vacuum fluctuations" ) and secondly to consider instead 
of the transition probabilities as function of E the sequence of moments 
of this function. When proceeding to arbitrarily short measurement times 
(while increasing the interaction coupling suitably) these moments will, in 
general, still diverge; however starting from the zeroth moment which stays 
finite one can subtract from the higher moments "perturbations" by the 
lower moments in such a way that the resulting (modified) moments all stay 
finite when sending the duration of measurement to zero. In this limit the 
modified moments are exactly what is measured by the balanced derivatives 
"in timelike direction" , i.e. by the balanced derivatives with £ tangential to 
the trajectory of the detector. 

Finally it will be shown that although these balanced derivatives do not 
span the whole iS^-spaces, their linear combinations can be used to obtain 
the relevant thermodynamic observables that so far appeared in the model 
of the massless, neutral Klein-Gordon field [3]. 

2 Specification of the model 

The quantum field considered is the massless neutral Klein-Gordon field on 
Minkowski spacetime (M 4 ,7y), described by a star-algebra A generated by a 
unit element 1 and elements </>(/) depending linearly on / £ Cq°(M 4 ) with 
the additional relations 



• <Kf)* = 4>U) 



• <£(□/) = o 



. W)M\=E(f,g)\ 



where 




2 see section 5 below for the notion of being at rest relative to a thermal state 
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and G a dv and G re t are the advanced and retarded fundamental solutions of 
the Klein-Gordon equation, respectively. 

States u of the field are given by functionals on A which in addition fulfill 
for all A e A and f 1} . . . f n G Cg°(R 4 ), 

• u{A*A) > 

• w(l) = 1 

• (/i,...,/ n ) i — > Lu((j)(fi) . . . 4>(f n )) is a distribution (depends continu- 
ously on the fi). 

There is a distinguished vacuum-state, denoted by Woo, and furthermore the 
states considered here are Gaussian Hadamard states, i.e. are determined by 
their two-point function C^(R 4 ) x C£°(R 4 ) B (f,g) uj{<p{f)(f){g)) which 
is such that (/, g) \— > uj{(f){f)(t){g)) — Ld 00 ((j)(f)(j)(g)) is a regular distribution. 
The physical picture of the measurements considered here is the following: 
An ensemble of quantum-mechanical detectors (two-level systems) moves 
along a (common, classical) trajectory 7, parametrized by proper time r, 
with each member initially in its ground-states. At some time the detectors 
are (smoothly) switched on, interact with the field for some time and are 
then switched off again. Finally the number of detectors in the excited state 
is determined, which yields the transition probability for a single detector 
(which is of course the same as the expectation value of the transition rate 
times the interaction duration). 

Mathematically, the free two-level detector system is described in the Heisen- 
berg picture by a two-dimensional complex Hilbert-Space (TId,< v >d) 
spanned by the two orthonormal states ip g and Ve of the detector. These 
two states are assumed to be eigenstates with eigenvalues zero and e of 
the detector-Hamiltonian Ho and furthermore the existence of a (time- 
dependent) self-adjoint operator r 1— ► M{t) := e lTHo MQe~ lrHD such that 
|< ip e , MqiP 9 >| 7^ is assumed. 

For a given state u of the quantum field, the coupled detector- field system 
is described in the interaction picture in the Hilbert-space TL^ ®TLd where 
TTwi Qw) denotes the GNS representation [T2l Chap. 4] of A belonging 
to uj. The initial state $ of the coupled system is taken to be tt u (g> ip g and 
the time evolution of this state is determined by 

id T $(r) = H int {T)$(r) := [x(t)0( 7 (t)) ® M(r)) $(r) 

where 7 : R — > M 4 is the detector-trajectory as described above, \ £ <5(R) 
is the detector switching function (real-valued and normalized by the re- 
quirement f R x( T )dT = 1) and R 4 3 x 1— > <j>{x) is related to ^(^(Z)) by 
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7Toj (</>(/)) = / (j>(x)f(x)dx in the sense of quadratic forms on (a subset of) 
H u . 

As one is interested in the transition probabilities of a (weakly) coupled 
detector without back-reaction effects a perturbative calculation seems ap- 
propriate; for a rigorous justification of this approach in a related situation 
see [6]. To first order perturbation-theory, the state of the detector-field 
system at time r is given by 

$( T ) = ^ ® V 9 - i / x(r') (0(7(t'))^) ® (M(r')V s ) dr' 

J — oo 

For the probability of finding the detector in the excited state Ve and the 
field in a state at large times, one then has 

|(* n <g>V> e ,$(oo)}| 2 = f [ { X (T')cl>( 1 (T'))n uJ ®M{T')lP g ,y n ®lP e ) X ... 

JR 

... x (*„ ® Ve, x(t")^(7(t"))^ ® M(r")^)dr / dr" 
and by summing over a complete set of ^ n in and using 

< ^ e ,M{r)^ g >=< ip e ,e iTH Moe- iTH ^ g >= e ltT < ^M^> a > 

the probability of finding the detector in the excited and the field in any 
state is 

P w (e) = |< ip e ,M ^ g >| 2 / / xVM^e-W- 7 ") x . . . 

JR JR 

... x < 0( 7 (r'))^, <A(7(r"))^ > dr'dr" 

= m [ [ X (r')x(r")e- ie(T '- T " ) a;^(7(T / ))^(7(T")))dT'dr" 
Jr Jr 

2 

where the constant of proportionality m := \< ip e ,Moip g >| depends on 
details of the detector but not on the field configuration. 
Now instead of comparing directly the transition probabilities -FL(e) and 
P L0/3 (e) in two states, one can in principle compare their difference to the 
transition probability P Wre / i n a common reference state. Choosing u)^ as 
iv re f this amounts heuristically to "removing the vacuum fluctuations" and 
the resulting (difference in) transition probability is then 

Pr(e) = ml I e- ies F u ,(T, S )x(r + s/2) X (T-s/2)dsdr (1) 

JR JR 

F u (j,a) := w(^( 7 (r + a/2))0(7(r-s/2)))-... 

...-^oo (<f>(rf (r + S /2))0( 7 (r-s/2))) 
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Assuming uj to be a Hadamard state, the integrand is smooth and com- 
pactly supported and therefore e i— > P^ en (e) is a rapidly decreasing, smooth 
function. Thus all moments of this function are defined and we turn to the 
analysis of these moments in order to obtain a means for local investigation 
of states. 

3 Convolution and moments 

As already mentioned above, in the case of global equilibrium and for x 
approaching a constant function, the dependence of the transition probabil- 
ities on e approaches a function that gives information about the thermal 
properties of the state under consideration. Disregarding for a moment the 
r-integration in (pQ), more rapid switching of the detector can be seen to 
disturb this function by convolution with a function that becomes wide as 
X becomes narrow, as is of course to be expected due to time-energy uncer- 
tainty. There is however a way to get around this, if one is only interested 
in the moments of this function and knows the moments of X- 
To see this, consider two rapidly decaying functions f,h £ S(M.). Denote 
the k-th moment of a function / E by 



and assume .M°[/i] = 1 (i.e. the integral over h is one). Then by a direct 
computation one has for the convolution f * h: 



For the special case k = this relation gives M°[f * h] = M[f], so the 
zeroth moment of / is always known once the zeroth moment of f * h is 
known. Noting that on the rhs of equation ([2]) the k-th moment of / is 
multiplied by one and the further terms involve only lower moments of /, 
one can therefore determine the moments of / from those of / * h and h by 
the recursion 





3=0 




M k [f] = M k [f *h]-J2(. )M j [f]M k - j [h] 



3=0 




(3) 



starting from A^°[/]. 
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4 Modified moments of transition rates and ele- 
ments of S x 

Now returning to ([I]) , for general x the application of the idea of the preced- 
ing section to the transition-probability is complicated by the r-integration. 
This problem can however be overcome by choosing x to be a Gaussian of 
width a centered at To: 

1 (t-tq) 2 



aV2ir 



Then Xa,r (r + s/2) X a,r (r - s/2) factorizes into X^/^/^X^o^) and 
the moments of e i— ► P^ en {e) are: 



M k [Pr] = m / / F w (r lfl ) X<T/ vs tT6 (T)dTe- to x^ ir ,oW^dc 




Now identifying the Fourier transform of s h- > \/2vr J K F u (t, ■5)x o .y v ^(r)dr 

with / and the Fourier transform of V^^X^o- with /i, this is the situation of 
the preceeding section and by the recursive procedure described there, one 
can obtain the fc-th moments of / which will be called V k . Explicitly they 
are given by 

Vt = mf [ F^s)x a/ ^ T0 (r)dre-^dse k de 
JM JM. 

By partial integration using the decay-properties of Xa,T and the regularity 
of F w this can be expressed as 

Vt = mf (-id s ) k FUT,s)\ s=oX<7/V 2 >7X> {r)dT 

JM. 

For this expression one can proceed to the limit a — ► 0, which shows that the 
V k are objects which can be measured over arbitrarily short time- intervals. 
In the limit one has 

limp* = m(-id s ) fc i^(T ,s)| s= o 

a— >0 

= m (-i/2) fe 5 s fc ( W -a; oo )(0(7(ro + S ))0(7(ro-5)))| s= o (4) 

Defining for a (sufficiently differentiable) function / : I 4 x I 4 -> I and 

ui, . . . , u n £ R 4 the balanced derivative g( Ml >->«n) j( x ) by 

n n 

Q(^-,u n ) f{x)=dti ^ dt j^ x + Y^t iUi , x-Y,t j u J )\ tl= ... =tn=0 (5) 

i=i j=i 
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one can rewrite ([!]) for inertial detectors with 7(70) = x and 7(70) = u as 
limP* = m(-i/2) fc g( M ' -' M )u;(x) 

w(x,y) := uj{(j){x)(j){y)) — 0J oo {4){x)(t){y)) in the sense of distrib. 

The definition © slightly generalizes the definition for gw* 1 !— >A*n) from [1] 
which is recovered by choosing Uj = ?7 ww e w with eo, . . . , e3 the basis vectors 
of M 4 , so the expectation values of elements m(—i/2) k g(°>— >°) : 2 : from 
the S^-spaces in the state uj used there can be interpreted as measurements 
of V k in the limit of arbitrarily short interaction of detector and field. A 
prominent example from this class is the Wick-square itself, which gives the 
expectation-value of the local temperature squared. 

5 Moving detectors and comparison with the full 
<S x -space 

As established in the last section, the measurement of balanced derivatives 
: can be described by a limiting process involving measure- 
ments carried out on an ensemble of detectors moving through the spacetime 
point x with a four-velocity u. These balanced derivatives however only gen- 
erate a subset of the 5 x -spaces and important (local thermal) observables 
like the (thermal) stress-energy tensor are not among them. 
Now global equilibrium states on Minkowski spacetime are not invariant 
under the full Poincare group but single out a set of inertial frames that 
only differ by rotations and translations and physically correspond to the 
observers being at rest with respect to the "gas" described by the thermal 
state. As thermality properties like the principle of detailed balancing hold 
only for those systems coupled to the field which are at rest in these inertial 
frames, for the investigation of an equilibrium state with an unknown as- 
sociated rest frame one should use not one detector, but detectors with all 
possible velocities smaller than the velocity of light relative to a given one. 
The detector behaving according to the principle of detailed balancing then 
indicates the rest frame of the given state, and starting from this one can 
then check whether the readings of the other detectors are compatible with 
the interpretation of being in relative motion to a thermal state. 
Whereas for a global equilibrium state, whose rest-frame is usually known a 
priori this discussion might sound rather odd, the hydrodynamical descrip- 
tion of gases by a velocity field varying in space and time can be rephrased 
as the statement that at each point the state looks like a thermal state 
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with reference frames at different points being in relative motion to each 
other. As this dependence of the frames on space-time is not known a priori 
but rather one of the informations an LTE- formalism should yield, it seems 
therefore sensible to not just consider one detector with a worldline passing 
through a spacetime point x but the set of all detectors passing through it. 
By the above procedure, one then obtains as local thermal observables the 
balanced derivatives g( u >--->") ; eft 2 : (x) for all time unit vectors u. By multi- 
plying with scalars, the requirement of normalization for u can be dropped, 
and forming linear combinations one furthermore has 

Ig( e i. e i) ; (ft 2 : = g( e O + |ej 5 eo+5 e j) : (ft 2 : + g( e 0-| e j: e o-|e j ) . ^2 . _ _ 

. . . -g( e °- e o) : eft 2 : 
2g( e o. e j) : (f) 2 : = g( e o+5 e i> e o+l e j) . ^ 2 . _g( e o-5 e i> e o-| e j) . 2 . 

g( e j. e fe) : 2 . = g(eo+|e i +|e A! ,eo+ie ;( +ie fc ) . ^2 . _ _ 
_g(e +ie J -ie fc ,e +ie j -ie fc ) . ^2 . _ 
_g(eo+|e fc ,eo+|e fc ) . ^2 . + g(eo-|e fe ,e -|e fc ) . ^2 . 

where ^ j, k,j = 1, 2, 3 and eo, . . . , e$ are the basis- vectors of M 4 which 
shows that balanced derivatives g(^) : cj) 2 : for fi, v = 0, ... ,4 can be ex- 
pressed as linear combinations of balanced derivatives g( u '") : (ft 2 : with 
different, timelike u (so that e.g. the expectation values of the thermal 
stress-energy tensor at x can be determined as a linear combination of mea- 
surement results of detectors moving through x with different velocities). 
In the models considered so far it was also possible to locally determine 
the entropy current and the particle density by introducing an approxi- 
mation process in the 5 x -spaces adapted to the notion of local thermality 
used. For a detailed description and motivation see [4]; at a technical level 
this approximation criterion boils down to the following: Given a function 
V + z) j3 i — ^ 6 K describing the (known) dependence of a macroscopic 
observable H on the (parameters (5 labeling) global equilibrium states, try 
to approximate S in the seminorms 

sup|S03)|=:||S|| B (6) 
pes 

(B a bounded set contained inside the open forward lightcone V + ) by the 
"thermal functions" of the elements in S x . The thermal function of an ele- 
ment from S x is defined as the mapping that associates to j3 the expectation 
value of this element in the global-equilibrium state belonging to (3; those 
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associated to 9^ : (ft 2 : (x) will be denoted by M . 

Investigating the global equilibrium states for the neutral, massless Klein- 
Gordon field (labeled by a timelike four- vector (3 already used above), the 
functions belonging to the particle density and the entropy current are found 
to be, respectively, a solution of the wave-equation and a gradient of such a 
solution (wrt. the parameter (3) 0j, [3]. That such functions can indeed be 
approximated in the above sense is established in [3]. The approximation 
given there only needs thermal functions l^Q^ corresponding to 9">— >') : (ft 2 : 
with I lightlikqfl, whereas here an approximation with thermal functions 
UfiQ^ of elements Q( u >—> u ) ; (ft 2 : with u timelike is desired. 
Relying on these results, it is thus sufficient to show that the functions l^Q^ 
for given, lightlike I can be approximated in the seminorms ([6]) by functions 
u (J 0' 1 with timelike vectors u. 

For the massless, neutral Klein-Gordon field, the thermal functions u^Q^ 
of elements B^ u >—> u ) : <ft 2 \ are proportional to [4] 

for m := even, and zero otherwise. This m-fold directional derivative 
is given by a sum starting with the term m\2 m (uf3) m ((3 2 )~ m ~ 1 and further 
terms proportional to (u 2 ) k (uj3) 1 (/3 2 )~ n (k, n £ N \ {0},/ G N) that vanish 
for u lightlike. Now for given lightlike I choose u = l + Seo. The first term in 
the sum is then itself a sum of m\2 m {lf3) m {f3 2 )~ m ~ l and functions bounded 
on B times powers of 5 and the remaining summands also give functions 
bounded on B multiplied by powers of 8 greater than zero (because of the 
terms (u 2 ) k , k G N \ {0}). For given B,l and e one can therefore choose 5 
such that (I + Se^^O** approximates / /X /J ' in the || • ||#-norm up to an error 
e. 

Summarizing, for the model of the massless, neutral Klein-Gordon field the 
relevant local thermodynamic observables can all be obtained by linear com- 
binations and limiting processes involving the balanced derivatives from sec- 
tion four, whose operational significance has been established. 
As far as the case of a massive Klein-Gordon field is concerned, the results 
up to and including the determination of the thermal stress-energy tensor 
using detectors in relative motion do not depend on the assumption of a 
vanishing mass; the question whether one can still approximate local ob- 
servables for the particle density and the entropy current in the case of a 

3 In this section fi is always assumed to be a multi-index and the Einstein summation 
convention is used 
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massive field is however more difficult to answer, as in this situation one 
has on the one hand to do the approximation in a different topology (for a 
definition and reasons for its use see [9]) and on the other hand the thermal 
functions of the generators of S x are more complicated. 

6 Conclusion 

Starting from an Unruh-de Witt detector model, it was shown that some 
modified moments characterizing the dependence of the transition rate of 
this detector on the separation of its levels survive the limit of going to 
arbitrarily short detector-field interactions (rescaling the strength of the 
interaction in the natural way). The value of these moments obtained in 
the limit coincides for geodesic detectors in Minkowski spacetime with the 
expectation values of the balanced derivatives defined in the context of non- 
equilibrium thermodynamics. Like other observables defining properties of 
quantum fields at a point (for example the stress-energy tensor), the objects 
have to be carefully defined in order to remain meaningful in the point- 
limit. Here this definition involves 

• Looking at differences between transition probabilities in different states 
instead of looking at the transition probabilities itself. 

• "Decoupling" time and frequency domain by choosing Gaussian switch- 
ing functions, choosing moments of the transition probability and re- 
moving the perturbing effects of lower on higher moments by the re- 
cursion-relation ([3|). 

The last point bears some similarity to a renormalization-group procedure: 
When going to higher moments, the high-frequency (small scale) properties 
of the field are emphasized and at the same time the low-frequency part is 
discarded in order to be able to proceed to a pointlike limit. 
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